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Abstract 

We find a new family of non-supersymmetric numerical solutions of IIB supergrav- 
ity which are dual to the Af = 1 cascading "conifold" theory perturbed by certain 
combinations of relevant single trace and marginal double trace operators with non 
infinitesimal couplings. The SUSY is broken but the resulting ground states, and 
their gravity duals, remain stable, at least perturbatively. Despite the complicated 
field theory dynamics the gravity solutions have a simple structure. They feature the 
Ricci-flat non-Kahler metric on the deformed conifold and the imaginary self-dual 
three-form flux accompanied by a constant dilaton. 



1 Introduction 



The problem of finding a holographic description for QCD is an outstanding chal- 
lenge. Even in the strongly coupled regime, when the dual background is expected to 
be weakly curved and potentially can be described by supergravity, to find the cor- 
responding gravity background is beyond our reach. Partially this can be attributed 
to the lack of supersymmetry. So far supersymmetry was the main vehicle to find 
new solutions and to guarantee their stability. Because of supersymmetry the dual 
solutions had a simple and elegant structure, but what we could learn about the dual 
field theories was bounded to supersymmetric dynamics. 

In this paper we want to make a step in the direction of finding a gravity back- 
ground dual to a non-supersymmetric field theory. As we are aiming at a confining 
theory the deformed conifold geometry of the Klebanov-Strassler (KS) background [1] 
seems a natural starting point. The KS geometry admits a global SU{2) x SU (2) sym- 
metry, which, although redundant from the QCD point of view, drastically simplifies 
the story on the gravity side. We would have to maintain this symmetry throughout 
the paper. 

At the next step we would like to break supersymmetry by means of perturbing 
theory by relevant operators, say, by giving mass to gauginos. This will break su- 
persymmetry while preserve the global SU(2) x SU(2). The corresponding gravity 
solution is straightforward to find when the gaugino mass is infinitesimally small and 
SUSY is softly broken [21 El Hj. Conceptually to find the corresponding solutions 
when the masses are large is also straightforward. Because the global SU (2) x SU (2) 
symmetry is preserved the corresponding lOd IIB supergravity background can be 
described by the so-called Papadopoulos-Tseytlin (PT) ansatz [5] that involves only 
ten functions of radial variable. The supergravity equations of motion reduce to ten 
coupled non-linear second order ODEs which one would need to solve imposing proper 
boundary conditions in the IR and the UV. 

In practice such a solution can be found only numerically but even in this case it 
is a formidable taskj^] To make the problem more manageable one needs to employ 
some trick which would drastically reduce the complexity of the system. For instance 
this was successfully done in [7] where a new family of gravity solutions dual to the 
baryonic branch of field theory was found. Thanks to supersymmetry the problem 

1 This question is successfully attacked in [5J. 
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reduced just to two coupled first order ODEs with simple boundary conditions while 
other eight functions from the PT ansatz are algebraically dependent on the first 
two. Later these solutions were generalized by use of a solution generating technique 
based on string theory dualities [HI El HO]- The resulting solutions are still attractively 
simple but develop various singularities. Typically these solutions have singular met- 
ric near the tip. And whenever solution is IR regular there is a deviation from the 
asymptotically AdS behavior in UV signaling the presence of the irrelevant dimen- 
sion eight operator. These solutions are of no use for us as we aim to find a regular 
non-supersymmetric background. Moreover we want to preserve the decoupling limit 
and therefore we do not allow perturbations by irrelevant operators. 

A good idea would be to preserve the simplicity of the supersymmetric solution 
while breaking supersymmetry explicitly. One example of such a trick in the context 
of D7-brane embedded into a supersymmetric background was developed in [11]. It 
was shown there that in a special case of constant dilaton background the supersym- 
metry conditions for D7 can be partially relaxed. Then the supersymmetry will be 
completely broken but the resulting solution for D7 will maintain the simplicity of a 
fully supersymmetric one. It is interesting to note that such a trick fails if the D7- 
brane is embedded into a more complicated supersymmetric background with running 
dilaton [T2]. In our case the situation is similar. Our trick is based on the following 
observation [T3l El EE]. The simplest version of the bulk supersymmetry condition 
requires the RR four form to be related to the warp-factor C4 = /i~ 1 Vol(IR 3 ' 1 ) and 
the three-form flux to be a (2, 1) primitive form. Furthermore, in the absence of 
D-branes the dilaton is constant while the unwarped 6d metric is Kahler and Ricci- 
flat. A remarkable observation is that if the three-form flux is an arbitrary imaginary 
self-dual (ISD) form (not necessarily (2, 1) primitive) this is enough for the fluxes to 
decouple from the metric/dilaton equations. Without sources the dilaton will have 
to be constant and one would end up with a Ricci-flat unwarped 6c? metric which is 
not necessarily Kahler. A non-Kahler metric and a generic ISD flux are both break- 
ing supersymmetry. We will denote this type of solution with the ISD three-form 
flux, a constant dilaton and a Ricci-flat Qd metric as the GKP (for Giddings, Kachru, 
Polchinski |15j ) background. These solutions are significantly simpler than the generic 
ones as the fluxes completely decouple from the equations for metric. 

Relaxing the SUSY condition only for the three-form flux while keeping the metric 
intact in the context of the KS background does not yield any new non-singular 
solution [2]. Therefore we are motivated to find a new Ricci-flat, presumably non- 
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Kahler, metric on the deformed conifold. Below we argue that the conventional 
Kahler Ricci-flat metric on the deformed conifod [16] - the unwarped metric of the 
KS solution - can be generalized to a one-parametric family of IR-regular Ricci-flat 
metrics which asymptote to the conventional metric in the UV. These new non- 
Kahler metrics can be used to construct novel gravity backgrounds which are regular 
everywhere and approach the KS solution in the UV. Therefore these backgrounds 
must be dual to the original N = 1 cascading "conifold" theory perturbed by some 
particular combinations of SUSY-breaking relevant and/or marginal operators. We 
do not expect these particular combinations of couplings to be special in any way from 
the field theory point of view. Quite the opposite, there must be a larger space of 
SUSY-breaking relevant deformations which lead to well-defined theories with stable 
vacua. 

To illustrate the point it is convenient to parameterize the family of new back- 
grounds which we find below, by U - the vev of the bottom component of the 
£7(l)baryon multiplet. This is because the supersymmetric solution for infinitesimally 
small U is of the GKP type [TTj . Hence U is a good coordinate to parametrize the 
new one-dimensional family of the GKP backgrounds, at least locally near U = 0. 
Value of U will specify the location of the new non-SUSY vacuum on the baryonic 
branch of the original SUSY theory. Let us start with the original SUSY theory in 
a vacuum on the baryonic branch with some non-trivial value of U. For finite U 
the dual gravity background is quite complicated: it is a Si7(3)-structure solution 
with running dilaton etc. [7]. Then we start turning on SUSY-breaking relevant and 
marginal couplings such that the vacuum value of U stays intact. Overall there is a 
large space of such combinations of couplings and the corresponding theories which 
we denote by M. v . Since SUSY is broken the corresponding gravity backgrounds 
presumably are even more complicated than the original point of our journey - the 
SU(3) structure solution of [7j with the given U. Yet, we claim, there is at least one 
particular point in JAu such that the gravity dual admits the simple GKP structure - 
it has an ISD flux, a constant dilaton and a Ricci-flat metric. These are the theories, 
and their gravity duals, we study in this paper. The schematic picture is shown in 
Fig. [TJ 

This paper is organized as follows. In the next section we numerically construct a 
one dimensional family of gravity backgrounds that corresponds to the blue line in fig. 
[TJ We construct the novel SU(2) x SU (2) invariant Ricci-flat non-Kahler metrics on 
the deformed conifold and then dress them up by the appropriate fluxes and warping. 
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Figure 1: We define M.y as a space of combinations of relevant /marginal couplings - 
and the corresponding theories - such that the bottom component of the {/(l)baryon 
multiplet has the vacuum value U. Schematically we denote x to be the coordinates 
on My (we do not specify the dimension of M.u) such that x = corresponds to 
the unperturbed M = 1 theory. Red dashed line corresponds to the motion in the 
space of theories starting from the unperturbed Af = 1 theory such that vev U stays 
constant. Most of these theories have complicated gravity dual. At a specific point 
the gravity dual accidently becomes simpler - it is of GKP type. The one-dimensional 
family of such theories forms the thin blue line. These are the theories we study. 

Then in Section [3] we discuss stability of the resulting solutions and their meaning 
from the dual field theory point of view. We outline our results in Section |4j 

2 Supergravity solutions 

2.1 Ricci-flat metric on the deformed conifold 

The starting point of our journey is the observation that the Gubser-Herzog-Klebanov 
(GHK) solution that corresponds to the vacua of the original M = 1 theory with the 
infinitesimally small U has the GKP structure [I7J ITS] : 

• Constant dilaton, e* = g s ; 

• Imaginary self-dual (ISD) three-form flux, iG 3 = *6G 3 , where * 6 is the 6d Hodge 
dual and G 3 = F 3 + -H 3 ; 
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• RR four-form C^ = h 1 Vol(IR 3 ' 1 ), where h is the warp factor: 

£j 2 dsj = h- l/2 dx^ + h l/2 dsl h ; (2.1) 

• Ricci-flat 6c? unwarped metric on M 6 . 

Beyond linear order the baryonic branch solutions of [7] do not preserve these prop- 
erties. A natural question then is wether it is possible to continue the GHK solution 
beyond linear order such that the Ricci-flatness of the unwarped metric and other 
nice GKP properties of the solutions are preserved. We will argue that this is in- 
deed possible and there is a one-dimensional family of such backgrounds which form 
the blue line in Fig [T] touching the horizontal line, the baryonic branch of the SUSY 
theory, at the origin. 

More broadly we can pose a question of finding non-singular Ricci-flat deforma- 
tions of the conventional Kahler metric on the deformed conifold. As an extra con- 
dition we will require the new metric to approach the original Kahler metric in the 
UV. For simplicity we also preserve the 577(2) x SU(2) symmetry. The most general 
metric compatible with these symmetries is given by the PT ansatz and depends on 
four functions of the radial variable r 



e^dsl = \e-^(dr 2 + gt) + e^\ C oMy){e\el + el) + e-\e\ + el)) 

-2sinh(y)(e 1 e 1 + e 2 e 2 ) V (2.2) 



We refer the reader to [2] for the definitions of the angular forms ej,6j, and g$. For 
the KS solution one finds z(t) =0 and: 



e Wp(r) = R 3^ sinh ( r ) ^ e 6g(r) = ^^4/5^) sinh 8/5( r ) ? ^(r) = tanh H_ 

8 V 2 



r „ , (sinh(2r) - 2r) 1/3 , . 

Wtae K ^= 2V3sinh(r) • (2 - 3) 

This Qd metric is a solution of the 10g? supergravity EOMs from [5], together with 
zero three-form fluxes and $ = const, provided we redefine some of the metric fields 

2 It is possible to introduce another fifth function A(t) by substituting tjj —> ip + A while keeping 
dtp intact. Obviously the constant part of A is a pure gauge associated with the action of U(1)r. 
Introducing A will not lead to any new solutions as the Ricci-flatness requires A = 0. 
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in the following manner (the expression for Ap^ ensures that the warp factor is one, 
h(r) = Ifl 

(e*,eV 2A ) pT = (2- 1 / 2 3- 1 / 4 e 2p+39 ,3 1 / 4 e p -",2- 1 / 2 3- 3 / 4 e 59 ) 1 



here 



e z 

e 9 = — — -, a = e z tanh(y) . (2.4) 
cosh(y) 

There is a Z 2 symmetry X that exchanges the two two-spheres of the conifold (and 
also changes sign of the three-forms). In the coordinates z\ satisfying £\ z\ = e 2 it 
acts by — >■ —24. This symmetry leaves the variables p(r), ?(r) and y(r) invariant 
and flips the sign of z(t). In the z = case the geometry is X invariant as in [1] (see 
also |2], HSJ). 



For the metric (2.2) to be Ricci-flat the functions p(r), q(r), y(r) and z(r) have to 
satisfy four second order ODEs and a first order constraint (the zero energy condition). 
Hence the full space of SU(2) x SU(2) invariant Ricci-flat metrics on conifold is seven- 
dimensional. Obviously the shift r — > r + const is a symmetry which we can fix by 
choosing the tip of the geometry to be at r = 0. Another simple parameter is the 
overall re-scaling of the metric given by q — > g+const. Alternatively, this shift controls 
the size of the three-sphere at the tip of the conifold which can be identified with the 
deformation parameter e. Thus we are left with only five non-trivial parameters and 
e. 

The counting above agrees with the linearized analysis in the vicinity of the sin- 
gular conifold z i = 0- The warped product of singular conifold and the Minkowski 
space is a gravity dual to a particular CFT [20] ■ Therefore all small Ricci-flat pertur- 
bations of the unwarped metric are in one-to-one correspondence with couplings and 
vevs of certain operators in field theory. There are three such operators in this case 



1. The gaugino bilinear A1A1 — A 2 A 



2^2 j 



2. The bottom component of {W 2 W 2 \ 



+ 1 

21 ID2I 



3. The bottom component of the ?7(l)baryon multiplet U ~ Tr(|A | — \B |) . 

We will denote these operators by their dimensions: 3, 6 and 2 respectively. The first 
two operators are even, while the last one is odd under X. Three operators correspond 



3 Our z(t) coincides with the one used in [T7] only at the linear order. Also, for z = our 
notations reduce to those of [2]. This is the reason for using p in the new notations. 
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to six small Ricci-flat perturbations of metric on the singular conifold. This is not in 
contradiction with our previous finding because in the singular conifold case the shift 
of t = ln(r 3 /e 2 ) and an overall rescaling of conic metric dr 2 + r 2 ds^ 1A coincide. Hence 
we are left with six independent linear modes that correspond to the three operators 
in the dual field theory. 

The analysis of small Ricci-flat perturbations around the deformed conifold should 
give similar results simply because in the UV region the deformed conifold approaches 
the singular one, up to 1/r corrections. Indeed there are six infinitesimal Ricci flat 
perturbations, two X odd and four X even. The X odd perturbations correspond to 
the coupling and the vev of U. The four X even perturbations correspond to the 
couplings and vevs of the operators of dimension 6 and 3. Because there are two 
different dimension three operators in this theory - two combinations of gauginos 
AjAi ± A2A2 - which mix in the deformed conifold case (but not in the singular 
conifold case) we will not be able to explicitly identify the operator responsible for 
this metric perturbation. In what follows we will simply denote this operator as AA 
and call the corresponding coupling the gaugino mass m\\. Four X even modes plus 
two X odd modes give six in total, but this is not in contradiction with five non-trivial 
parameters we found above. This is because one of the six modes is actually the shift 
of t i.e. the r-derivative of the background solution. 

These six infinitesimal Ricci-flat perturbations around the conventional conifold 
metric have been found explicitly. Unfortunately the full description is quite tech- 
nical. We briefly mention it below but a detailed knowledge of these modes is not 
necessary to follow our logic. An uninterested reader can skip until the next para- 
graph. The X-odd Ricci-flat perturbations were found by GHK in [T7]. The X-even 
ones can be found using the formalism to study linear X even SU(2) x SU(2) in- 
variant perturbations around the KS solution proposed by Borokhov and Gubser [21] 
and later developed in [25J. The full space of solutions is parameterized by sixteen 
constants X\, . . . , X 8 , Yi, . . . ,Y 8 . To find the Ricci-flat perturbations of the conifold 
metric one can impose that full background is of the GKP type, i.e. the dilaton is 
constant, three-form flux is ISD etc. Then ignoring an overall rescaling one recov- 
ers exactly four modes. In the notations of [4j these modes are parameterized by 
X 3 , X4,Y 2 ,Y 3 while X 2 = —2X 3 /3, Y\ = —5Y 2 /3 and all other X^ and Fj are zero. 
Only three modes are truly non-trivial as the mode Y 2 = const, Y 3 = X 3 = X 4 = is 
the derivative of the background with respect to r. It can be attributed to the vev 
of an operator of dimension 3 in a sense that the linear mode asymptotes to 1/r 3 at 
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infinity. Then X 4 corresponds to the coupling m\\ while X 3 ,Y 3 correspond to the 
coupling and the vev of the dimension six operator correspondingly. 

Since we know the linear Ricci-flat perturbations explicitly we can analyze their 
behavior in the UV and IR. There is only one UV divergent mode which corresponds 
to turning on the coupling of the irrelevant dimension six operator in the dual field 
theory. To keep the UV asymptotic the same as in the conventional case and hence the 
field theory well defined in the UV we want to keep this mode "turned off" . In the IR 
there are two regular modes. One is X odd and corresponds to the infinitesimal motion 
along the baryonic branch. The other is the X even mode that turns on the dimension 
six coupling - the same mode we discussed above. So there is no Ricci-flat IR and UV 
regular X invariant infinitesimal deformation of the conventional deformed conifold 
metric. To find a new Ricci-flat background one has to consider the X breaking modes 
as well. 

The two IR regular linearized modes can be extended into two-parameter family of 
IR regular Ricci-flat metrics on the conifold. With vanishing fluxes, constant dilaton 
and h = 1 the IIB supergravity action of [5] reduces to the 6c? Hilbert-Einstein 
action, and so we can read from it the four Ricci flatness equations. We relegated 
these equations to Appendix A. The two-parameter family of the IR regular solutions 
is giver^jby: 



r+ I -- + — <2 |r 3 + ... 



3 V 45 
1 15 / 7 45 



q{r) 4 V 8 26 

,M = ^ + (-I. 6 -l.<f-g.« 1 )^ + .... (2.5) 
We write here only the result for q(r) since adding constant to q(r) corresponds to 



the overall rescaling of the metric (2.2). This rescaling can be in turn absorbed in 
a redefinition of the deformation parameter e. Fixing of e in the new solution is an 
important issue and we will specifically address it in the end of this section. 



4 We choose the field parameterization here (e 10p ( T ) instead of p(r), q 1 (t) instead of q(r) etc) so 
that all the expressions will be regular both at the origin and for large r. 
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Figure 2: The red points show the values of (i and (2 for which the Ricci-flat metric 
is regular in both the IR and the UV. Notice that the points near the origin lie on the 
parabola (the blue line). The deviation of the points with larger d from the parabola 



most likely indicates the contribution of the (f term in (2.7). 



We introduced here two parameters (1 and (2 such that the X symmetry flips the 
sign of (1 while keeping (2 invariant. The point (Ci 5 C2) = (0,0) is the conventional 
deformed conifold metric, i.e. the unwarped metric of the KS solution, while the 
infinitesimal (2 and £i correspond to "turning on" the dimension six mode and the 
GHK mode respectively. Notice that both £1 and (2 appear in the sub-leading order 
in p(r), q(r) and y(r) and so the solution is indeed regular at r = for any Q2 and 
(1 exactly as the original KS background. 

Let us stress again that we solve the full-nonlinear Ricci-flatness equations and 



(2, Ci are n °t assumed to be infinitesimally small in (2.5). 



Regularity in the UV provides only one constraint on (j 2 and therefore we expect 
to find a one-parameter family of non-singular Ricci-flat metrics on the deformed 
conifold, which approach the conventional metric at large radius. The boundary 
conditions at r = 00 should coincide with the deformed conifold values (2.3): 



,W P (t) 



1 9 

W) W 2' 



1 



0. 



(2.6) 
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Figure 3: The four lines show numerical solutions for e 10p ( T \ l/q(r), e y ^ and z(r) for 
(d,C 2 ) = (0.193(9), 0.009) (the last point on the previous plot). 

In Appendix A we study the subleading terms and show that the modes of dimension 
2, 3 and 6 indeed appear in the large r expansion as expected. 



Near the origin the leading perturbation (see the z(r) expansion in (2.3)) is pro- 
portional to d (meaning it is I odd) and therefore at least for small £2 and £i is it 
convenient to parameterize the family of new Ricci-flat metrics by (\. 

Ca(Ci) ~ Ci + 0(Ct) , (2.7) 

where the term is ruled out because of the X-parity. We found the function £2 = 
C2(Ci) numerically using the shooting technique and confirm a generic prediction that 



the expansion (2.7) starts with £f (see Figure [2J). We also present a numerical solution 
for the functions e 10p<T \ l/g(r), e y<yT \ z(t) for some particular values of £i and £ 2 , see 
Figure [3] Finally, on Figure [4] we compare our numerical results for e y ^ with its 
analytical conventional (KS) counterpart of |16j . 

It would be interesting to know how far one can move along this new family of 
solutions and if there is any limit on the resulting value of U. These questions would 
require careful numerical studies which we leave for the future. 
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Figure 4: The blue solid line show the numerical result for e y ^ with (Ci 5 C2) — 
(0.193(9), 0.009), while the red dashed line corresponds to e yo ^ = tanh(|) of the 
conventional deformed conifold of [16J. Notice that red line converges to the asymp- 
totic faster than the blue line. This is so because for large r the function 1 — e y ( T ' 
goes like e~ T for the the deformed conifold, and as e _T / 3 for our solution. 



Ricci-fiat metric on deformed conifold: a summary 

The construction of the novel family of Ricci-fiat metrics on deformed conifold pre- 
sented above constitutes a mathematical result which could be of interest in its own 
right, with no connection to the ten dimensional supergravity or dual gauge theory. 
To make this aspect of our work more accessible to reader less familiar with the holo- 



graphic context we briefly outline here main steps of section |2.1| Our goal was to 
construct novel Ricci-fiat metrics on the deformed conifold - the complex three-fold 
described by the equation Ylt z f = ^ ■ We construct the Ricci-fiat metric by help 



of the ansatz (2.2). This ansatz is by no means general. In particular it explicitly 



preserves the SU (2) x SU (2) symmetry. Since all functions p, q, y, z depend only 
on the radial coordinate r the Ricci-flatness condition = reduces to a set of 



coupled ordinary differential equations (A.l). Our ansatz (2.2) is invariant under the 



shift t — > t + const. We choose this freedom to set r = as the origin - the tip of 
he conifold. That would be reflected in the appropriate behavior of the metric there 
i.e. particular boundary conditions for p, q, y, z at r = 0. We find a two-parameter 
family of solutions regular at r = 0, i.e. satisfying these boundary conditions, by ex- 
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panding in power series in r for small r. The result is given by (2.5). Two integration 
constants £1,2 are free parameters at this point. We want the metric to approach the 
usual conical form at large radius r — » 00. Here by usual we mean the canonical form 
of the metric on the cone over T 1 ' . There is one Ricci-flat metric on the deformed 
conifold, found by Candelas and de la Ossa [TB], which does just that. This particular 
metric is captured by our ansatz. At the next step we study the behavior of linear 
fluctuations around this metric, all within our ansatz of choice. The system of linear 
equations reveals that there is only one mode which tends to grow at large r. All 
other modes actually vanish in this limit. On general grounds we expect that by 
choosing appropriate infinitesimal (1,(2 one can fine-tune the behavior of the metric 
near r — > such that the growing linear mode will not be present at large r. Hence 
the behavior at large r would be the same as in the solution of [16], as well as the 
linear analysis around it: there is still only one growing mode at large r. In this way 
we find that there must be a line in the (1 , (2 plane such that for each point on this 
line the resulting metric has the same asymptotic behavior at large r. These metrics 
are Ricci-flat by construction, as well as regular everywhere including appropriate 
behavior at the tip r = and at infinity r — > 00. 



2.2 Adding fluxes and warping 

In the previous subsection we outlined the way to find a regular Ricci-flat metric 
on the deformed conifold such that it approaches the conventional metric at large 
radius. The next step is to supplement such a metric with a smooth ISD three-form 
flux which would approach the KS asymptotic in the UV. Below we give an argument 
why this is possible for any metric found in the previous section. 

We start with the PT ansatz for the three-form flux which automatically includes 
M units of RR flux through the three sphere at the tip and write down the ISD 
condition: 

/ = (cosh(y) cosh(z) + sinh(?/)) 2 — F • ( cosh(2?/) + smh(2y) cosh(z)) 
k = cosh 2 (?/) sinh 2 (,2) + F ■ ( cosh(2y) — sinh(2y) cosh(z)) 

F = \{k~f) (2.8) 

X = — cosh(y) sinh(2;) ( (cosh(j/) cosh(z) + sinh(?/)) — 2F ■ sinh(y) ] . 



13 



These functions are related to the functions used in [S] in the following way: 



(h h h 2 , b, X ) PT = {-Pg,{k + f), -Pg s (k - /), 2F - 1, 2P x ) here , (2.9) 

where P is a constant proportional to M, P = —M£ 2 S /A. It is easy to verify that for 
e y(r) _ tanh(r/2) and z(r) = this system of equations reduces to the one discussed 
originally in pQ. 

Naively the solutions are parameterized by four integration constants. One of the 
constants is a shift of \ — > \ + const which corresponds to a gauge transformation 
B — >■ B + const ■ dg 5 and hence is unphysical. Another constant is a shift (/, k) — > 
(/ + b,k + b) for a constant b which corresponds to a "large gauge transformation" 
B — > B + b ■ u>2, where w 2 is a Betti form on the base of the conifold. This shift 
introduces extra Page D3-brane charge to the system. This symmetry will be fixed 
by the requirement that the Page D3-charge vanishes at the origin, namely that 
-BA-F3 = 0atr = 0. Therefore there are only two non-trivial parameters left. 

In the UV the metric and fluxes should approach the KS values. Linear analysis 
around the KS background shows that apart from the simultaneous shift of / and 
k there are indeed two modes [2]: one is UV regular and the other is UV singular 
diverging as e T . The former is related to the A = 3 operator, while the latter is the 
(0,3) flux dual to the A = 7 operator / d 2 d(W 2 W 2 ) + . Eventually we would like to 
"turn off" this irrelevant operator "killing" one of the two free parameters. 

It seems like we get a one-parameter family of UV-regular solutions of the ISD 



equations (2.8). This is not the end of the story, however, as we have yet to address 



the IR regularity. Differentiating /(t) and k(r) we can arrive at a second order PDE 



for F(t) only. Using (2.5) we see that, exactly as in the KS case, near r = the two 
solutions of this equation are r _1 and r 2 . Plugging the latter into the equations for 
/(r) and k(r) one finds that these functions have the same regular behavior in the 
IR like in p]. This in turn implies that \F 2 \ = \H 2 \ does not diverge at r = and 
so the background is IR regular. To be more specific, requiring that /(r) and k{j) 



vanish at r = we derive from (2.8) and (2.5) that there is a one dimensional family 



of IR regular ISD fluxes parametrized by 

^(t) « (1 + Gi) ~ , *(t)«(1 + C S )-|, X^-^Ci(4 + C 3 )-r 2 . 

(2.10) 

In particular £i, ( 3 = for the KS background. 
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Figure 5: The solid blue and the dashed red lines correspond to the numeric solution 
for F(t) with (i = 0.193(9) and the same function in the KS geometry, respectively. 



We now have to find £3 so that the divergent A = 7 mode is "turned off" in the 
UV. Then for large r we have F(t) ~ 1/2, while both f(r) and k(r) go like r/2, 
reproducing the UV behavior of the KS solution. Employing the shooting technique 
once again we found that £3 = 0.186(8) for (1 = 0.193(9) (the last point in Figure 
[2]). The corresponding numeric solution for F(t) is shown in Figure [5j With F(r) at 
hand, we can easily calculate the remaining three functions. In Figure [6] we provide 
a plot for x{ T )- As a matter of convenience, we fixed the integration constant so that 
x(t) vanishes for large r, see Figure [6j 

To complete our discussion of the ten dimensional supergravity background we 
need to accompany our solution by the warp factor of the 10c? metric (2.1). The 
equation for h(r): 

fo poo 

h(r) = ^fa / dxe~ 4p{x) - 6q{x) ((l - F(x))f(x) + F(x)k(x)) , (2.11) 

where: 

a = 4(g s Mf s ) 2 ■ e" 8 / 3 . (2.12) 
This equation immediately reveals that h(r) is regulai|^]at the tip and approaches the 
5 We fixed the integration constant in the equation for h{r) by requiring that the warp function 
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Figure 6: The numerical solution for x( T ) with Ci = 0.193(9). 



KS behavior at large r. 

The last touch is to fix the deformation parameter e (as we have already mentioned 



below (2.5) this is equivalent to adding a constant to the field q(r)). 

Since e is the only dimensional parameter of the solution we can only talk about 
fixing it in the context of comparing e for two different solutions. We want to fix e(£i) 
such that the corresponding family of gravity backgrounds would correspond to the 
same dual field theory perturbed by some relevant operators in the UV. In particular, 
this would require running of the Maxwell D3-charga^] Q ^3 with the scale to be the 
same at large radius r for all £i- A convenient reference point is £i = 0, i.e. the KS 
solution. Thus, to match the KS theory in the UV we require: 

QD 3 (r c ) = QZ§{r c ) (2.13) 

at a sufficiently large cutoff radius r c . Here the Maxwell D3-charge is (see, for example, 



vanishes at infinity exactly like in the unperturbed KS background. 

6 To be more precise, one has to compare the running of the Page charges, since unlike the Maxwell 
charge it is quantized and so has the proper interpretation as the gravity dual of the cascading gauge 
theory rank. It turns out, however, that the matching of the Maxwell charges yields exactly the 
same result for e. 
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my- 



Q m {r) = ^M 2 ((l - F(t))/(t) + F{r)k{r)) . 



71 



(2.14) 



It seems that this expression is completely independent of e. Notice, however, that 
t = ln(r 3 /e 2 ) and so requiring that the two Maxwell charges in (2.13) coincide for 
the same r = r c (but different r's) we can read from (2.13) the ratio between the KS 
parameter exs an d the parameter of the perturbed theory, e. For instance, for the 
theory corresponding to the last point on Figure [2] (the one with d = 0.193(9)) we 
have Q D3 {t = 11.866(2)) = Q$§(t = 12) which implies e/e KS = 1.069. It is worth 
emphasizing here that for large enough r c (or equivalently r) the value of e/exs is 
not sensitive to the variation of the cut-off, since asymptotically Qd3( t ) ~ r + const. 



A more detailed discussion of the condition (2.13) and the resulting constraint on 
e in the context of SUSY and non-SUSY vacua of the KS theory can be found in [I] . 



3 Dual field theory and stability 

In this section we discuss the field theories dual to our gravity backgrounds. We 
suggest that these non-supersymmetric theories correspond to the KS theory per- 
turbed by both a single and a double-trace operators. We also argue that, at least for 
large M, these theories are long-lived both perturbatively and non-perturbatively. 



3.1 The dual field theory interpretation 

Once we outlined the way to find the gravity backgrounds we want to analyze their 
meaning from the dual field theory point of view. As we require the backgrounds to 
approach the KS solution in the UV the resulting field theory must be the original 
cascading N — 1 "conifold" theory perturbed by some relevant and/or marginal 
operators. In the single trace sector these are the X-odd bottom component of the 
[/(l)baryon multiplet U and the gaugino bilinears. It is easy to see that the marginal 
top components of W± are not turned on: the value of dilaton and running of the 
Maxwell D3-charge are the same for large r as in the KS case. (In any way adding 
top components of W± would only result in the renormalization of gauge couplings). 
Hence the solutions in question are dual to the vacua described by the Af = 1 KS 
theory perturbed by some, not necessary small, combination of the gaugino masses 
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7Ti\ 1 \ 1 ± Tfi\ 2 \ 2 and the operator U. We will not be able to identify the particular 
combination of mj^j ± m\ 2 \ 2 that is getting turned on (see our discussion in the 
previous section) and in what follows will refer to this combinations simply as m\\. 
The vacuum solutions also acquire a non-trivial vev of U, as well as vevs of other 
relevant and irrelevant operators. 

For most operators it is straightforward to distinguish the vev from the coupling 
as the two modes have different powers of 1/r where the radial coordinate 

r = e 2 /V/ 3 . (3.1) 

The operator U, however, has dimension two and therefore the two modes mix. The 
corresponding wave-function, which for large radius can be defined through: 



sds 2 = \/W)^l + 4-^1- 4)*u (3.2) 

has the following asymptotic [TT] : 

*u = /?ln f^j +a . (3.3) 

Here ro is some IR scale which we choose such that, in the absence of multi-trace 
deformations, the coefficient (3 corresponds to the vev of U while a corresponds to 
its coupling. In general ro may depend on various couplings of the theory, such as 
ttl\\. In the KS case, when all relevant couplings (including m\\), are turned off, 

ro = e 2/3 e l/3 [J7]_ 

We do not know r (m\\), although this can be established in principle. Therefore 
even if we find our solutions numerically with a very high precision, we will not know 
a(Ci)F] Nevertheless, on general grounds we expect: 

a ~ d 3 + 0(Ci 5 ) , P ~ Ci + 0((1) . (3.4) 

This scaling is easy to explain. f3 corresponds to the vev of U which, by definition 
of Cij scales as £i- Moreover at the linear order in C\ the solution coincides with the 
GHK background, i.e. has zero coupling a. Since a is odd under X symmetry, as well 



7 In practice the main difficulty here is to find the solutions numerically with enough precision 
to distinguish small a which appears only in the £i order from the large logarithmic contribution 
PMr/ro) ~ Cir. 
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as £1, it can only start with a cubic term. Let us demonstrate that the same scaling 
is compatible with (3.3). From the fact that m\\ is Z-even we find that: 



m AA ~ Ci + 0(Ci 4 ) • (3-5) 

Now, we generally expect r to be mAA-dependent such that dr /dm\\ ^ 0. The 
change in r$ can be reabsorbed into a which indeed indicates that a ~ (f . 

A non-zero value of a implies that besides the gaugino mass the dual field theory 
is perturbed by the dimension two operator U . Definitely this can not be the full 
story: the theory perturbed by the potential V = all does not have a vacuum at 
small U. We find a consistent interpretation only at the next order in U - at least for 
small d the gravity solutions in question can be interpreted as a dual to the theory 
perturbed by the marginal double-trace^] operator clI 2 /2. The coupling c will alter 
the boundary conditions in the UV and hence the dual field theory interpretation of 
a. More precisely, when the field theory perturbed by a potential: 

U 2 

V = aU + c—, (3.6) 



the boundary condition for tyjj becomes [261 

a = a + c(3 . (3.7) 

When c = we get the usual relation a = a, i.e. in this case a is the coupling of 
U . Since the theory with the V = all potential has no vacuum state there should 
be no IR regular solution with a = a ^ 0. In the general c ^ case the resulting 
field theory has a unique vacuum at U = —a/c. The corresponding gravity dual 



background should have a and (3 such that they satisfy (3.7) and the geometry is 
regular in the IR. Since the IR regularity requires a = 0, we see that = —a/c, in 
full agreement with (3 = U and dV/dU = 0. 

It is important to note that the gravity solution itself can not distinguish between 
different combinations of a, c which result in the same vev of UP\ In fact the same 
gravity solution admits interpretation as both being dual to a field theory with and 
without the double trace deformation (see [281 [29] an d [30] correspondingly). Using 
this freedom we would like to interpret our solutions as the gravity dual of the field 
theory perturbed by certain single trace relevant operators and also the marginal 



8 The dimension of the double trace operator U 2 is twice the dimension of U up to 1/N corrections. 
9 We thank I. Klebanov and J. Maldacena for discussing this point. 
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double-trace operator U 2 . Our logic here was to show that it is a consistent interpre- 
tation when in the dual field theory a ~ Ci an d c ~ Ci • 

Strictly speaking, the consideration above and the equivalence between vacuum 
value of U in the field theory and in the bulk is only valid as long as the theory with 
a = c = has a flat direction for VEV of U. This is indeed true for infinitesimal d 
as it it can be shown (using the charges under R-symmetry) that the gaugino masses 
can only contribute to the potential for U at the subleading order m\ x ~ (f. But 
we believe that our interpretation - that the solutions we found are dual to the KS 
theory perturbed by gaugino masses and the potential V = aU + cU 2 /2 - should be 
correct beyond the infinitesimal order in (i. In this case the gaugino masses induce a 
potential for U (on top of V — aU + clI 2 /2 which we add "by hands" in the UV) and 
the calculation of vacuum value of U is obscured on the field theory side. Similarly 
on the supergravity side (3.7) still holds but we do not know which combination of 
a, (3 corresponds to the IR regular background. Thus we can not easily compare field 
theory and gravity to prove our point. At the same time the fact that we included 
all relevant /marginal operators (allowed by symmetries) into consideration and the 
fact that this interpretation works for small (i, suggests that the same interpretation 
should be correct beyond the infinitesimal order in (i. 



3.2 Stability 

The field theory interpretation suggested above ensures the perturbative stability for 



infinitesimally small Ci as long as c > (see (3.6)). There are many field theories 
dual to the same gravity background and we choose the one with positive c: for any 
theory with a, c dual to a given background the theory with —a, — c is dual to the 
same background as well. Indeed when a = this flipping will not modify (3, as one 



can see from (3.7), and so the gravity solution will be the same. 

On the gravity side the stability would follow from the analysis of small pertur- 
bations around the classical solutions in question with the boundary conditions: 

5a = c 5/3 . (3.8) 

Clearly the sign of the perturbation mass squared will be sensitive to the sign of 
c. More specifically, in solving the equations for massive modes, the sign of the 4d 
dimensional mass squared will be sensitive to the relative sign of 5a and 5(3, which 
in turn is fixed by c. 
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By interpreting the gravity solutions as being dual to the field theory perturbed by 
the double-trace operator we ensured perturbative stability, at least for the solutions 
with infinitesimally small (j. For larger (\ we would need to calculate the mass-squre 
mfj of small perturbations of U to check that is it positive. For infinitesimally small 
(i, mfj ~ c and the next correction appears in the (f order 

ml ~ d 2 + 0((t) • (3.9) 

Since the coefficient in front of the first term is positive, mfj is positive at least for some 
finite range of £i < 1. Hence the vacuum states, and the dual gravity backgrounds, 
are perturbatively stable for £1 < 1. 

Besides the perturbative stability, one may also worry about possible tunneling 
between the 2M vacua of the original KS theory with different phases of A 3 . When 
(i is small, i.e. when m\\ is small, one can estimate the shift of vacuum energy in 
each of those vacua as [2]: 

V ~ MRe (m\\A 3 ) . (3.10) 

Our solutions correspond to m\\ aligned with A 3 - they have the opposite phase 
such that their product is real. This is so since we have not included in our setup 
the mode A(r) mentioned briefly in footnote [2] and its counterpart in the fluxes. It 
remains unclear, though, whether ?7?,aaA 3 is positive or not. Hence our solutions may 
be unstable non-perturbatively. 

Let us show that in any case the theories in question are long-lived, if we do not 
depart too far from the KS solution. At the leading order the tension of the domain 
wall separating various vacua remains the same as in the supersymmetric case [31 J : 

T n = M\A\ 3 \l-e inn/M \ , n = l,...,2M. (3.11) 

When m\\ is small enough we can use the thin wall approximation and estimate the 
size of a new vacuum bubble by minimizing the following effective action: 

S huhhlc ~ AVR 4 - T An R 3 . (3.12) 

For the adjacent vacua with An = 1 we have T ~ |A| 3 while AV is 1/M suppressed. 
For a one leap transition into the true vacuum An = M we get T ~ MA 3 and 
AV ~ Mm aaA 3 . The decay rate for tunneling into a new vacuum with A 3 — > A 3 e J7r / M 
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is smaller than the decay rate for tunneling into the true vacuum in one leap. The 
latter is suppressed by: 

^bubble _ , 

e -M(A/m AA ) 3 _ (3.13) 

Because of large M for m AA < A i.e. £i < 1 the solutions in question are metastable. 

This analysis of stability may be too conservative. In fact it may be that the signs 
of Tfi\\ and A 3 are aligned such that the solutions we find correspond to the true 
vacuum. 



4 Outline 



We constructed a one- dimensional family of IIB supergravity solutions which are 
dual to the M = 1 gauge theory perturbed by some combination of relevant single 
trace and marginal double-trace operators. Thus the family of solutions corresponds 
to a line in the space of couplings <?i(Ci)- The origin of the line ^ (0) = corresponds 
to the original M = 1 theory pQ . Except for the small region near the origin d — > 
the couplings are not infinitesimally small and the corresponding operators explicitly 
break supersymmetry. Hence our solutions are gravity duals for non-supersymmetric 
field theories. 

The original M = 1 theory has only one massless mode which acquires a positive 
mass, at least for the theories within a range d < 1. Therefore the solutions in 
question are perturbatively stable. There is a possible non-perturbative instability 
associated with tunneling into a vacuum with A 3 ew = —A 3 . But at least for d < 1 
the decay rate is exponentially suppressed. 

The geometry of the solutions is extremely simple. The unwarped metric on the 
deformed conifold is Ricci-flat but not Kahler for d ^ 0. The three-form fluxes are 
ISD and the dilaton is constant. There is also a warp-factor which has a similar 
behavior to its KS counterpart. The solutions are regular and all fields approach 
their KS values in UV, up to 1/r corrections. 

The apparent simplicity of the solutions makes them a natural arena to study 
strongly coupled dynamics of the non-SUSY confining gauge theories. Similarly to 
the KS background, which was a playground for numerous phenomenological models, 
our solutions can be used for model-building with an advantage that the resulting 
models are not supersymmetric. Because of the GKP structure of the solutions the 
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probe D3-brane in such backgrounds experiences no force. Hence our solutions can be 
a natural starting place for the models of stringy inflation based on the dynamics of a 
mobile D3-brane [321 |33j [3U |35j [361 122] . In a similar way the GKP structure makes it 
possible to use the trick of [11] in its original form and embed the U-shaped D7-branes 
inside the conifold. Eventually one arrives at a holographic model of baryonic matter 
with the SUSY broken explicitly by the relevant operators in the gauge sector, and 
not only by the flavor sector. Moreover because of the KS asymptotic at infinity the 
resulting baryons will develop a realistic attractive potential at least in some range of 
parameters [37] ■ Finally, it would be interesting to understand the appearance of the 
flat moduli associated with the mobile D3-branes from the point of view of the dual 
non-supersymmetric field theory (see [38] for a similar situation in a 3c? dimensional 
theory) . 
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Appendix A A. The Ricci flatness equations 



Here we summarize the Ricci-flatness equations for the 6c? metric (2.2) 

2 



p + p{Ap + 6g) 



15 



e -2o P (i + 2 cosh 2 y sinh 2 z) 



(A.l) 



1 



q + q(Ap + 6g) 



H e 10p cosh y cosh z H — sinh 

15 5 

Q 

e -2o P (i + 2 cosh 2 ?/ sinh 2 z) 
2 



135 
16 



H e 10p cosh u cosh 2 sinh z y 

45 15 



y + 2/(4p + 6g) 



cosh y sinh y ( -e 20p sinh 2 z - -1 



4 _ 10p cosh2 



coshy 



+ (l + i 2 )) 



and: 



z + i(4p + 6q + 2y tanh y) = sinh z ( —e 20p cosh 



z - te- 10p — 



(A.2) 



3 coshy / 

These equations can be derived from the SU(2) x SU(2) action of [3] with vanishing 



fluxes and the relations (2.4): 



S 



dre 4 P +6 9 



15p 2 



^ 2 + l {f + cosh 2 (y)z 2 ) 



(A.3) 



+ Qe" 20p (l + 2 cosh 2 (?/) sinh 2 (^)) - 2e" 10p cosh(y) cosh(^) + ^ sinh 2 (?/)^ 



The solutions of these EOM are also subject to the zero energy (zero Hamiltonian) 
condition following from the action. This condition can be used to find q(r), and so 
we are left only with the equations for p(r), y{r) and z(r). 

To study the large r behavior we have to expand the equations linearly around 



(2.4). The equations for Sy(r) and 6z(r) are both homogenous. The 5z(r) equation 



is solved by e~ 2r / 3 and re~ 2r,/3 , as expected for the the A = 2 operator. The two 
solutions for Sy(r) are e~ T and e~ T / 3 corresponding to A = 3. Finally, the homogenous 
part of the equation for Sp(r) gives two modes, e~ 2r and e 2T//3 , implying that the mode 
has dimension A = 6. 
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